Light-matter interaction in Mie nanolasers 
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By deriving a three dimensional vector set of Maxwell-Bloch equations, we report on an ab-initio 
investigation of a spherical Mie nanolaser. Parallel numerical simulations predict a rich physical 
scenario, ranging from a nontrivial vectorial energy-matter interplay in the initial stage of evolution 
to mode competition and dynamical frequency pulling effects. Application of these effects could 
favor the realization of nonlinearly-controlled largely- tunable nanolaser devices. 
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Since the pioneering investigations of Gustav Mie in 
the early days of the last century [1], the theoretical and 
experimental study of spherical resonators has deserved 
much attention from the scientific community. This inter- 
est stems from the analysis of fundamental processes such 
as scattering, energy propagation through disordered me- 
dia and cavity quantum electrodynamics, and from the 
large number of applications in photonics, chemistry, me- 
teorology, astronomy and sensing [3, 0, 0, 0, El, 0, S, • 
Mie resonators display ultrahigh quality factors Q 
10^) [l^, and are able to confine and store electro- 
magnetic energy for long times within small dielectric 
volumes; Mie whispering galleries modes (WGM) can 
sustain enhanced stimulated emission processes so that 
'thresholdless' lasing can be observed De- 
spite the widespread scientific production, the theoreti- 
cal understanding of Mie lasers is still an open field of 
research. For what concerns laser emission, existing the- 
oretical approaches either rely on a semi-classical treat- 
ment based on light interaction with a two level system 
[H, 0, [itI or on rate equations ^ , with numerical anal- 
ysis limited to one and two-dimensional media fio']. The 
former approach, as was recently observed [20] , is not rig- 
orous as the number of physical dimensions grows above 
one while the latter, accounting only for the contribution 
of atomic populations, misses important informations on 
quantum coherence [2l|, and cannot be effectively pur- 
sued to perform any realistic ah initio computation at 
ultra- fast time-scales (<ps). Furthermore, although low 
dimensional models permit a simplified analysis with re- 
spect to the fully vectorial set of Maxwell equations they 
left the general picture unknown. Which is the outcome 
of a strongly nonlinear ultra- fast and multi- dimensional 
interaction on scales comparable with the wavelength? 
Beside being central for the development of laser theory, 
this question has fundamental importance in quantum 
mechanics [22|, soliton [23] and chaos theory ^]. How- 
ever, to the best of our knowledge, neither ah initio inves- 
tigations were reported nor quantum-mechanical models 
of energy-matter interaction for three-dimensional nanos- 
tructured lasing systems were developed. Owing to the 



large interest in the study of nano-scaled particle aggre- 
gates [25, 2II, 15^, the theoretical investigation of such 
dynamics is of considerable practical importance. 
In this Letter we develop an ab initio rigorous theoret- 
ical model of light interaction in the presence of ampli- 
fying (or dissipative) materials, deriving a three dimen- 
sional vector set of Maxwell-Bloch (MB) equations within 
the real representation in terms of the SU{n) algebra 
(2^ . We discretized the resulting equations on a Yee grid 
and numerically solved them within the Finite-Difference 
Time-Domain (FDTD) method [29]. The MB-FDTD ap- 
proach is then applied to investigate the ultrafast dy- 
namics of Mie nanoresonators. Specifically, we perform a 
series of numerical experiments by investigating the pro- 
cess of laser emission from a single nanosphere, covered 
by a layer of active material, for different pumping rates. 
This theoretical analysis yields the following scenario: 

i) In the transient before the steady state regime strong 
interactions between all components of the electromag- 
netic field and the atom take place. Such a vectorial 
dynamics is the key ingredient for the reaching of a sta- 
ble steady state. 

ii) When the lasing process settles up, modes possessing 
frequencies close to Mie resonances get excited and com- 
petition dynamics begins. 

iii) As the pumping rate gets higher, a strongly multi- 
mode regime occurs. In this situation modes with higher 
quality factors tend to grow faster than the others and 
dynamical frequency pulling effects are observed. In 
particular, as the pumping rate increases the principal 
frequency of laser emission moves toward smaller wave- 
lengths encompassing large shifts; a process radically dif- 
ferent from the standard nonlinear frequency shift ob- 
served in lasers. The mechanism investigated here can be 
at the basis of the development of nonlinearly-controlled 
largely- tunable nanolaser devices. 

We begin by writing Maxwell equations in an isotropic 
medium: 
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being P the material polarization. The latter can be 
decomposed as the sum of a linear contribution Pl = 
eo(e^ — 1)E plus a nonlinear term Pnl = —^^aiQ) mod- 
eling the quantum-mechanical interaction between the 
electromagnetic field and the atoms being e the electric 
charge, Na the density of polarizable atoms and (Q) the 
expectation of the displacement operator Q with respect 
to the quantum state |'0(r, t)). To derive a quantum me- 
chanical model of the material, we consider a four-level 
system, with three degenerate levels, in which three elec- 
tric mutually orthogonal dipole transitions can be excited 
by linearly polarized electromagnetic waves with energy 
equal to the atomic resonance hcoo (Fig. [TJ-left). The 
total Hamiltonian operator H of the system can be then 
expressed as H = Hq + H/, with the unperturbed Hamil- 
tonian Ho of components Hij = huJoSij{Si^2 + ^i,3 + 
{6ij is the Kronecker delta) and the self-adjoint interac- 
tion Hamiltonian H/ = eE • Q. The displacement opera- 
tor is Q = G + G''' with Gij = qo^i,i{^j,2^-\-^j,3y-\-^j,4z)^ 
being qo the typical atomic length scale and x, y, z unit 
vectors along x, y and z axes, respectively. The dynam- 
ical evolution of the atomic system is expressed by the 
Liouville equation of motion of the density matrix op- 
erator p = \ijj){ip\: ih^ = [H,p]. By exploiting the 
symmetry of p under the SU (4) group [28] , we introduce 
the coherence vector S of components Sj: Sj = Trips j] 
{j = 1,2, ...,15) being Tr the trace operator and Sj the 
j— th generator of the SU{4:) algebra [30]. The time evo- 
lution of the coherence S is then: 

— =fS-7(S-S(°)), r,^ = -Tv{U[s„SH]), (2) 

being 7 a phenomenologically added [22[ diagonal matrix 
of nonuniform relaxation rates with 7^^ = 1/Tj, S^^^ a 
vector containing the initial populations of the system, 
r is skew and its nonzero components turn out to be: 
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with fti = [i G [x^y^z]) and p = ego- Finally, the ex- 
pectation value of the displacement operator Q expressed 
in terms of S is: 

(Q) = Tr [pQ] = go (x5i + y^4 + z^g) • (4) 

Equations Q-dl]) represent the three dimensional, full- 
wave vector set of Maxwell-Bloch FDTD equations mod- 
eling light-matter interaction in the presence of resonant 
quantum transitions. 
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FIG. 1: (Color online), (left) Energy diagram of a four- level 
atomic system with a triply-degenerate excited state; (right) 
sketch of a single Mie nanoresonator with a d— thick layer of 
amplifying material. 



Equations ©-([I]) have been discretized on faces and 
sides of a cubic Yee grid while ([Tj) on its center. Standard 
leapfrogging has been used as time marching algorithm 
for the FDTD, second order Crank-Nicholson for the so- 
lution of the MB set. A Uniaxial Perfectly Matched Layer 
(UPML) is employed to absorb outgoing waves. Follow- 
ing 0, we keep into account fluctuations by adding to 
the electric field a Markovian noise process with Gaus- 
sian statistic, the latter obtained from a Box-Miiller ran- 
dom number generator algorithm. The code has been 
parallelized within the Message-Passing-Interface (MPI) 
standard. 

We perform a series of numerical experiments by con- 
sidering a spherical resonator (Fig. [T]-right) of r=300 
nm radius, while taking = 8.41 as a realistic value 
for typical high-index materials (e.g. Ti02 ), covered 
by a layer of active material of thickness d = 100 nm, 
with the same refractive index of the interior part and 
transition at uoq = 3.88 • 10"^^ rad/s (Ao = 483 nm) in 
proximity of Mie resonances, as detailed below. By as- 
suming ^ 0.1 nm, we obtain a coupling coefficient 
p = 4.8 • 10~^^ C • m. The decaying constants were set 
to Ti = 10 fs {i = 1,2,4,5,6,7,9,10,11,12,13,14) and 
Tj = 1 ns {j = 3,8,15) according to the existing lit- 
erature [1^. The pumping rate is settled by the den- 
sity of polarizable atoms Na in the excited state, which 
is assumed to be constant in time and fixed by an ex- 
ternal continuous excitation involving additional energy 
levels, as for standard laser systems. Furthermore, ow- 
ing the small sphere radius, we assume a spatially uni- 
form pumping of the active layer. The used discretization 
Ax = Ay — Az = 10 nm and At — 0.006 fs guarantees 
both high accuracy (with 50 points per wavelength at 
A = 500 nm) as well as temporal stability. 
We begin our analysis at a low pumping rate [Na — 10^^ 
m~^). To investigate the transient regime, we collect at 
each time step all the electric field components E^^ Ey^ 
Ez^ sampled in proximity of the sphere center, and the 
output electromagnetic energy density f = E • D + H • B 
(Fig. [2]). Ai t ^ 600 fs the electromagnetic fields begin 
to coherently grow. Owing to the sphere-air high in- 
dex contrast, the electromagnetic energy is well confined 
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2: (Color online). MB-FDTD results for Na = 10^^ 
m (a)-(c) time evolution of the electric field components 
(sampled in proximity of the sphere center, to avoid symmetry 
effects), (d) energy density S (isosurface plot) with xy (down) 
xz (left) and yz (right) slices in the sphere middle plane. 
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FIG. 3: (Color online), (a) Time series spectrogram com- 
puted from MB-FDTD analysis for Na = 10^^ m"^ (the 
mode-competition dynamics is denoted by the beating mode- 
amplitudes); (b, solid line) Mie theory extinction factor Qext 
for a single nanosphere of diameter 600 nm and Cr — 8.41; (b, 
dashed line) active medium gain bandwidth. 



within the structure and the lasing mode experiences low 
losses (Fig. [2]). This stage of evolution is characterized 
by a mutual energy exchange between the atom and each 
electromagnetic component (Fig. [2^-c). The latter will 
eventually interact on the fs scale (fixed by Ti = lOfs) 
with each other owing to the presence of nonlinear non 
diagonal coupling terms in the density matrix ([3|) . From a 
physical perspective, energy exchange originates from the 
presence of a triply degenerate atomic excited state (Fig. 
[U left panel) , which yields a nontrivial vectorial interplay 
between the electromagnetic field and the atom. Such 
a contribution, ignored by previous theories, is the es- 
sential ingredient for an isotropic laser three-dimensional 
dynamics. Without a four- level atomic system, in fact, 
no vectorial interaction is possible and no realistic steady 
state can be theoretically predicted. 

Once the transient stage has been concluded {t ^ 2 
ps), multimode lasing settles up. Even if nano-sized, the 
sphere has high quality factors that do not prevent las- 
ing to occur (Fig. [3|). To study the time evolution of the 
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FIG. 4: (Color online). Frequency demodulation analysis: 
time evolution of modal amplitude ai(t) (a) and phases 0i,2(t) 
(c-d); (b) ID map obtained with a Poincare surface of section 
of (ai,a2) at ai = 180 kY/m. 



steady state regime we make use of a spectrogram (Fig. 
[3]) and frequency demultiplexing (Fig. [4]). The former 
shows a mode- competition dynamics between two modes 
in proximity of coq. In Fig. [3)3 we plot the extinction fac- 
tor Qext of the nanosphere as calculated from Mie theory 
which shows the available modes (Fig. [3)3 solid line) 
in the gain bandwidth (Fig. [3)3 dotted line). The two 
modes (as shown in Figs. [3^-b) correspond to the near- 
est WGM modes with respect to the medium resonance 
uq. To characterize the mode-competition dynamics, we 
extract amplitudes ai(t), a2{t) and phases (/>i(t), 02 (^) 
at carrier frequencies close to WGM modes by means of 
frequency demodulation (Fig. [4]). In principle, a four 
dimensional phase-space [ai, a2, (/>i, (/>2] supports a rich 
dynamics ranging from periodic motion to chaos. How- 
ever, at this low pumping rate, no synchronization 0] 
occurs and each phase independently linearly grows, ex- 
cept in a small set of time instants where is close to 
zero (Fig. [4^-d). This is the "free-run" regime, as origi- 
nally predicted by Lamb [l^, where the phase variables 
can be averaged out from the coupled-mode equations 
and only the amplitudes remain as dynamical variables. 
As a result, we are left with a two-dimensional phase 
space which, according to the Poincare-Bendixson theo- 
rem, can only sustain fixed points or periodic motion, as 
shown in Fig. [4^ and in Fig. [4)3, where we plot the map 
obtained from the Poincare surface of section of (ai, a2). 
The bi-modal steady state is characterized by a limit cy- 
cle where both amplitudes ai, a2 are trapped together. 

As the pumping rate Na gets higher, a large number 
of WGM is above the threshold for laser oscillation. To 
investigate the transition from two-mode to multi-mode 
regime, we calculate spectrograms obtained from MB- 
FDTD simulations for increasing pumping rates Na G 
j2o26 iq28] ^-3 (^p-g rpj^g dynamical evolution ob- 
served is characterized by two competing effects. As the 
frequency uj increases with respect to ujo- i) the gain 
experienced by the mode gets smaller; ii) the quality fac- 
tor Q of the WGM mode becomes higher. Therefore, 
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FIG. 5: (Color online). Spectrograms for increasing pumping 
rates: Na = 10^^ m"^ (a), Na = 10^'^ m"^; (c) Power density 
spectrum for various Na- 



although low wavelength modes experience small gain, 
they tend to grow faster with respect to modes at higher 
wavelengths. The result of this competition is a dynam- 
ical frequency pulling towards smaller wavelength as the 
pumping rate grows (Fig. [Sb). The predicted frequency 
shift is very large 100 nm, see Fig. [Sb) and quanti- 
tatively it depends on the morphology of the nanosphere 
resonance spectrum (see Fig. [Sb), which fix the involved 
wavelenghts and the corresponding Q-factors; in princi- 
ple it can be controlled by acting on the shape of the res- 
onator. Dynamical frequency pulling is then expected to 
play a fundamental role in the nonlinear dynamics of laser 
systems, either based on a single sphere or on ordered 
and disordered resonator ensembles. The induced large 
shift can be exploited to develop largely-tunable nano- 
sized laser devices operating at different wavelengths by 
acting on the input excitation; such applications would 
require a proper engineering of the cavity and will be the 
subject of future works. 

In conclusion, we have developed a rigorous theory of 
electromagnetic energy-matter interaction in the pres- 
ence of nonlinear resonant media. The derived equations 
are applied to investigate the process of laser emission 
of Mie nanospheres. By means of parallel simulations 
and frequency demultiplexing analysis, we discuss both 
the transient and the steady state regime of laser emis- 
sion. During the transient, a nontrivial mutual energy 
exchange between the electromagnetic wave and the res- 
onant medium occurs, thereby justifying a rigorous ab- 
initio investigation. On the other hand, the lasing state 
observed by increasing the pumping rate is characterized 
by a dynamical transition from two-mode to multi-mode 
regime, the latter accompanied by competition phenom- 
ena and dynamical frequency pulling. Such results have 
important implications in the theory of nano-lasers, in 
both ordered and disordered systems and are expected to 
stimulate new experiments and applications where single 



Mie resonators are employed, including nonlinear optics, 
colloidal physics, chemistry and quantum mechanics. 
We acknowledge support from the INFM-CINECA ini- 
tiative for parallel computing and S. Trillo for useful dis- 
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